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ABSTRACT 



Context. The origin of very energetic charged particles and the production of very high-energy (VHE) gamma-ray emission remains 
still a challenging issue in modern pulsar physics. 

Aims. By applying a toy model, we explore the acceleration of co-rotating charged particles close to the light surface in a plasma-rich 
pulsar magnetosphere and study their interactions with magnetic and photon fields under conditions appropriate for Crab-type pulsars. 
Methods. Centrifugal acceleration of particles in a monopol-like magnetic field geometry is analyzed and the efficiency constraints, 
imposed by corotation, inverse Compton interactions and curvature radiation reaction are determined. We derive expressions for the 
maximum particle energy and provide estimates for the corresponding high-energy curvature and inverse Compton power outputs. 
Results. It is shown that for Crab-like pulsars, electron Lorentz factor up to y ~ 10 7 can be achieved, allowing inverse Compton 
(Klein-Nishina) up-scattering of thermal photons to TeV energies with a maximum luminosity output of ~ 10 31 erg/s. Curvature 
radiation, on the other hand, will result in a strong GeV emission output of up to ~ (10 34 - 10 35 ) erg/s, quasi-exponentially decreasing 
towards higher energies for photon energies below ~ 50 GeV. 

Conclusions. Accordingly to the results presented only young pulsars are expected to be sites of detectable VHE y-ray emission. 
Key words. Pulsars - Acceleration of particles - emission processes 



1. Introduction 

One of the fundamental problems in pulsar physics is re- 
lated to the origin of the observed non-thermal emission. 
While it seems evident that efficient particle acceleration and 
emission processes must be operating in a pulsar's rotating 
magnetosphere, current theoretical approaches differ widely 
in their assumptions about the localization of the relevant 
zones: According to standard polar cap models, for example, 
charged particles are uprooted from the neutron star's surface 
by strong electrostatic fields (IRuderman & Sutherl and 1975l l. 
Still close to the star, these particles are then assumed to be 
efficiently accelerated along open field lines in parallel elec- 
tric fields induced by space-charge-limited flow (Mi chel 199 U . 
field-line curvature (lArons 1983) and/or inertial frame drag- 
ging effects (Muslimov & Tsygan 1992 1. In most cases, the par- 
allel electric field component is shorted out at some altitude 
by the onset of electron-positron pair cascades in strong mag- 
netic fields (one-photon pair production), either initiated by 



els, for example, a critical issue has always been the question 
whether particles can indeed gain enough energy inside the gap 
to account for the observed non-thermal radiation from y-ray 
pulsars. Several scenarios have been proposed to enlarge the 
gap zone and consequently increase the corresponding energy 
output (e.g.JUsov & Shabad 19851 lArons & Scharlemann 1 979 



Muslimov & Tsygan 1992| l, yet y-ray emission from Crab-like 



curvature! Daugherty & Harding 1982 ) or Inverse Compton ra- 
diation dDermer & Sturner 1994K Outer gap models, on the 
other hand, assume that primary particles are efficiently ac- 
celerated in vacuum gaps in the outer magnetosphere, induc- 



ing pair cascades through yy-pair production ( |Cheng etaL 1986 
Chiang & Romani 1994 Hirotani 20071 ) . In all these approaches 



the maximum attainable particle energy is either limited by the 
gap size or radiation reaction. In conventional polar cap mod- 
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pulsars still proves challenging to account for. The efficiency 
of particle acceleration along magnetic field lines has also been 
studied more recently based on numerical solutions of the struc- 
ture of a stationary, axisymmetric and force-free magnetosphere 
of an aligned pulsar (Contop oulos et al. 1999[ ). According to the 
results obtained, the relativistic magnetospheric outflow is not 
accelerated efficiently enough to account for the production of 
high energy gamma-rays. In some respects, this result may not 
come unexpected as the magnetic field configuration is restricted 
to be force-free, thus preventing efficient acceleration. In a rather 
different approach, Beskin & Rafikov (2000) analyzed the ac- 
celeration of a (stationary) two-component, electron-positron 
outflow in a monopole magnetic field configuration for high 
(Michel) magnetization parameters cr = B 2 l%jmm e c 2 , indicating 
the ability of the field to sling particles to high velocities. Here, 
n e is the electron density, B the induction of the magnetic field, 
c is the speed of light and e and m e are electron's charge and the 
rest mass respectively. Considering plasma dynamics close to 
the force-free regime (first-order correction), they showed that 
for small longitudinal currents very high Lorentz factors can be 
achieved, with almost all of the electromagnetic energy being 
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converted into the kinetic energy of particles (y max ~ <x) in a 
thin layer close to the light cylinder surface. 

In the present contribution we consider another accelera- 
tion mechanism that may help to overcome the energy prob- 
lem arising in some polar cap-type models. To this end, we 
explore the acceleration of co-rotating charged particles in 
an idealized monopole-like magnetic field region close to the 
light surface where the parallel electric field component is ef- 
fectively screened out, but where inertial (centrifugal) effects 
become important in describing the plasma dynamics. This 
follows earlier suggestions by Gold (1968, 1969) about ef- 
ficient particle acceleration close the light cylinder in a co- 
rotating neutron star magnetosphere (see also Ruder man 19721 1. 
A detailed analysis of centrifugal acceleration along rotating 
straight field lines in the test particle limit has been pre- 
sented by Machabeli & Rogava (1994), showing that due to 
the relativistic mass increment the radial acceleration of a 
particle changes sign, similar to results obtained for parti- 
cle motion close to a Schwarzschild black hole (Abramowicz 
& Prasanna 1990). Based on this, the plasma motion in 
pulsar magnetospheres has been analyzed and equations de- 
scribing the behavior of a co-rotating plasma stream have 
been derived (e.g., IChediaet al. 19961 IMachabeh et al. 20051 . 
More recently, the generalization to curved field lines (e.g., 
Archimedes spiral, where a particle may asymptotically reach 
the force-free regime) has been examined and the conse- 



quences of radiation reaction analyzed (Ro gava et al. 2003 
Dalakishv ili et al. 20 07 ). Independently, applications of cen- 
trifugal particle acceleration to milli-second pulsars were con- 
sidered and curvature radiation effects discussed in Gangadhara 
( 119961 ) (see also |Thomas & G angadhara 2007 for a recent gen- 
eralization). In a wider context, the efficiency of centrifugal 
particle acceleration was studied for Active Galactic Nuclei 
(AGN) ([Gangadhara & Lesch 1997| |Rieger & Mannheim 2000 



Osm anov et al. 20071 ), based on scenarios where AGN jets origi 
nate as centrifugally-driven outflows (Bland ford & Payne 1982[ ). 

In the present paper, we analyze the efficiency of centrifugal 
acceleration for Crab-like pulsars, taking constraints imposed by 
co-rotation, inverse Compton interactions and curvature radia- 
tion into account. The paper is arranged as follows: In Sec. 2 
the radial particle motion due to centrifugal acceleration effects 
is described and co-rotation constraints discussed. In Sec. 3 we 
examine possible radiative feedbacks on the process of accel- 
eration for typical millisecond pulsars, considering some major 
limiting processes: inverse Compton scattering, curvature radi- 
ation and pair creation. In Sec. 4 the relevance of our results is 
shortly discussed in the context of recent observational evidence. 



2. Centrifugal acceleration of particles 

2.1. A simplified approach 

We consider an idealized, single-particle approach, where in the 
local frame of reference, each particle is only subject to the ac- 
tion of the centrifugal force and gains energy while moving out- 
ward along the magnetic field. The field configuration is sup- 
posed to be almost straight, at least inside the co-rotation zone. 
This seems a reasonable assumption to make as efficient parti- 
cle acceleration is expected to take place on characteristic length 
scales much smaller than the light cylinder radius. When applied 
to open field lines, our results are thus not expected to be very 
sensitive to further magnetospheric details. If a denotes the an- 
gle between the magnetic field B and angular velocity Q, then 
the effective field line rotation is Q e = Osinc and the light 



surface is at a distance - r^o/ sin a, where r^n = cP/2n 
is the light cylinder radius and P = 2n/Q. is the pulsar pe- 
riod. The dynamics of a co-rotating particle can then be de- 
scribed by the following equation (Machabeli & Rogava 1994 
Rieg er & Mannheim 2000[ ) : 



d 2 r _ Q 2 e r 
dfl ~ 1 - QV2 
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and its Lorentz factor can be expressed as 
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where m - 1/y 2 , is essentially determined by the initial con- 
ditions. Hence, for a secondary pair plasma produced close to 
the neutron star with Lorentz factor yo ~ (10 3 - 10 4 ), for ex- 
ample, efficient centrifugal acceleration to high energies accel- 
eration can take place close to the light surface, i.e., in a layer 
Ar/rL ~ yo/y. Using Eq. (f2]), the characteristic timescale for 
centrifugal acceleration can be approximated by 



r 



dy/dt 4k sinarw'^y 1 / 2 



(3) 



Obviously, for a particle approaching the light surface, the ac- 
celeration timescale decreases with y and the Lorentz factor can 
increase dramatically unless co-rotation can no longer be main- 
tained or radiation reaction becomes important. 

2.2. Co-rotation constraints 

Suppose that the co-rotation zone extends outwards from the 
neutron star up to the vicinity of the the light surface (Gold 
1968, 1969). Because of strong synchrotron losses, electrons 
will quickly lose their relativistic perpendicular energy, i.e. on 
a timescale which for most pitch angles iff is much smaller than 
the transit timer, = y 2 sin 2 ip t s to their ground Landau state, 
so that they may be approximately described as moving one- 
dimensionally along the field lines. Yet, even if one neglects ra- 
diation reaction (e.g., curvature losses, see below) co-rotation 
will only be possible as long as the kinetic energy density of 
the electrons ynm e c 2 does not exceed the energy density in the 
field B 2 /8tt (Alfven corotation condition). For a number density 
n = M ncs, where M denotes the multiplicity (number of secon- 
daries to number of primaries) and «gj = QB cos a/ (2nec) 
0.07 Bcosa/P [particles cm 3 ] the classical Goldreich- Julian 
number density close to the star, the co-rotation condition im- 
plies an upper limit for achievable electron Lorentz factors of 



y^e^XlO 5 



BP 



Mcos a 



(4) 



For a Crab-type pulsar with P ~ 0.033 s and B 10 6 sin 3 a G (at 
the light surface r^, assuming an internal dipolar field structure), 
for example, this results in y^ x e 2x 10 1() tan a sin 2 a/M, while 
for a Is pulsar (B ~ 10 sin 3 a G at r\J) the upper limit would be 
of the order of y™^ e =s 4 x 10 7 tan a sin 2 a/M. These values sug- 
gest that for milli-second pulsars co-rotation may indeed last up 
to the very vicinity of the light cylinder. Note that close to the 
light surface, the number density required to screen out a po- 
tential parallel electric field component is much larger than «gj 
employed above, and in fact given by (Goldreich & Julian 1969) 



1 



'GJ,L 



«GJ 



(1 -Q 2 r 2 /c 2 )' 



(5) 
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The onset of a pair production front in polar cap models, close 
to the neutron star, is usually expected to result in a multiplic- 
ity M = n/n G} ~ 10 2 - 10 5 and a mean Lorentz factors of the 
secondaries of yo ~ 10 3 - 10 4 (e.g.,|Daugherty & Harding 1983 
IMelrose 19981 Baring 2004). On the other hand, using Eq. K2J 
Eq. (0 can be expressed as «gj,l — »gjT/To- Hence, in order to 
achieve electron Lorentz factors of, e.g., y ~ 10 7 via centrifugal 
acceleration, one requires M ~ y/yo ~ (10 3 - 10 4 ), which seems 
well possible given the parameter range above. Equation (0]i then 
suggests conversely, that for Crab-type pulsars, Lorentz factors 
up to y ma x,e ~ 10 7 can indeed be achieved if radiation reaction is 
negligible. Acceleration will then occur in a narrow region very 
close to the light surface with Ar/r^ ~ 1 /M, This seems remi- 
niscent of earlier results about particle acceleration in the pres- 
ence of small longitudinal currents (Beskin et al. 1983; Beskin 
& Rafikov 2000). 

3. Emission constraints 

In realistic astrophysical environments, radiation reaction will 
impose additional constraints on the efficiency of any particle ac- 
celeration process. For pulsars, important limitations could arise 
through inverse Compton scattering with ambient soft photons 
field or synchro-curvature losses along curved particle trajecto- 
ries. Pair production (i.e., one-photon or photon-photon) on the 
other hand, could possibly lead to a suppression of detectable 
high energy y-rays. 



3. 1 . Inverse Compton interactions 

Thermal radiation as well as synchrotron radiation by secondary 
electrons could in principle lead to a non-negligible target pho- 
ton field for Inverse Compton (IC) interactions and thereby limit 
achievable electron energies. 



3.1.1. Inverse Compton with thermal photons 

It has often been assumed that IC interactions with thermal 
photons from the neutron star surface are generally negligible 
far away from the surface because (i) the photon density de- 
creases with distance r and (ii) charges and photon are travel- 
ing in almost the same direction, so that (anisotropic) inverse 
Compton losses become exceedingly small (e.g., |Morini 198 it . 
While the first consideration is certainly true, the latter may 
not necessarily be the case. In fact, if electrons are co-rotating 
with the plasma, their main velocity component close to the 
light cylinder is expected to be in the azimuthal direction, im- 
plying a preferred interaction angle of almost 90 degree, so 
that IC interactions with thermal photon field may possibly be- 
come relevant for milli-second pulsars. Although pulsars are 
born at very high temperatures T ~ 10 11 K, their surface tem- 
peratures quickly cool down to ~ 5 x 10 6 K by various neu- 
trino emission processes and thermal emission of photons (e.g., 
ITsuruta et al. 2002t lYakbvlev & Pethick 2004) . Standard (mod- 
ified UCRA, plasma neutrino and photon cooling) models pre- 
dict a surface temperature above 10 6 K for pulsars with ages 
t = P/2P ~ 10 4 yr (neutrino cooling stage), and below 10 5 K 
for pulsars exceeding 10 6 7 yr (photon cooling stage). Using the 
standard cooling curve, one can employ an approximate phe- 
nomenological description for the temperature-age dependence 
given by (Zhang & Harding 2000} 



T(t) 



s /'10 6 yr 
5.9xl0 5 K' 



0.1 



forr < 10 52 yr. 



T(t) 2.8 x 10 5 K 



10 6 



yr 



0.5 



forr > 10 52 yr- 



(7) 



Hence, for a Crab-type pulsar (t 10 3 yr, rjr^fl — 10~ 2 ) (where 
r s is the stellar radius) the surface temperature may be approx- 
imated by T - 1.2 x 10 6 K. The Planck function then peaks 
at around v max = 2.8 IcT/h 7 x 10 16 Hz, corresponding to 
a photon energy e p h - 0.3 keV. IC scattering thus mainly oc- 
curs in the extreme Klein-Nishina regime. As a first order ap- 
proximation for the single particle (non-resonant) Klein-Nishina 
Compton power, one can employ the expression derived by 
Blu menthal & Gould (1970)] assuming a (quasi-isotropic) black 
-body photon distribution, i.e., 



(jj(m e ckT) 2 



c,KN - 



16fi 3 



ln W_ lA)81 )(^\ 



(8) 



where y is the electron Lorentz factor. This implies a character- 
istic IC cooling timescale close to of of 



ym e c 
he = — °c y , 

*c,KN 



(9) 



which to first order is proportional to y. Acceleration, on the 
other hand, occurs on a timescale f acc oc 1 /y 1 ^ 2 , so that for elec- 
tron Lorentz factors y ~ yo, IC cooling will not impose any con- 
straints on achievable particle energies. This suggests that for 
Crab-like pulsars electron Lorentz factors are essentially limited 
by co-rotation and not by IC radiation reaction (cf. eq.|4|i. If so, 
then detectable IC emission at ym e c 2 ~ 5 (y/10 7 ) TeV may oc- 
cur. We can roughly estimate the possible TeV luminosity by 
multiplying the corresponding particle number with the single 
IC power P c .kn, i.e., L^ v ~ ncjMAVP c kn> which gives 

L TeV „ 10 3l/_B_W T N 2 pq 

IC \10 5 G/\1.2x 10 6 K/ IsincJ s/ ' 

where AV = \nxr\d ~ 2nxr^/M, with ^ < 1 denoting the 
deviation from isotropy and d ~ (yo/2y) r L the thickness of 
the layer close to r^, in which the highest particle energies are 
achieved. The (pulsed, non-steady) TeV luminosity could thus 
be as high as ~ 10 31 erg/s, consistent with existing upper limits 
derived by current g r ound-based y-ray instr u ments (IACT) (e.g . , 
ILessard et al.~2000l lAharonian et al. 20071 lAlbert et al. 2 008 ), 
yet possibly accessible to the next-generation CTA-type instru- 
ments. 



3.1.2. Inverse Compton with infrared photons 

Secondary synchrotron emission could possibly lead to a non- 
negligible photon field in the infrared-optical regime where IC 
interactions may occur in the Thomson regime. For the Crab 
pulsar, the (isotropic, phase-averaged) near infrared-optical lu- 
minosity is of order L ~ 10 33 erg/s, turning significantly 
downward for lower frequencies (e.g., Middled itch et al. 1983 ; 
Eikenberry et al. 1997 ISollerman 2003b . This suggests a pho- 
ton energy density close to the light surface of order ~ 
L /(47itlc) ~ 10 5 sin 2 o' erg/cm 3 , comparable to the thermal 
one. Approximating the single particle (non-resonant, quasi- 
isotropic) Compton power by P c ~ cr T cy 2 u p h/(l + x), where 
x - ye p h/m e c 2 , the characteristic IC cooling timescale close to 
then becomes 



(6) fic - 



ym e c 



2 4.1xl0 7 sec 



r « P h,L 



:(l+2y/10 6 ). 



(ID 
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Comparing acceleration, occurring on f acc (Eq. |3j, with IC cool- 
ing, occurring on fic, implies 



y 



< 2x 10 



to 



(l+2y/10 6 ) 2 ro sin 2 a 



(12) 



and verifies that IC interactions with the infrared-optical photon 
field will not impose a severe constraint on the maximum achiev- 
able Lorentz factor. IC up-scattering by electrons with y ~ 10 7 
would again produce emission at around 5 TeV. If the observed 
NIR-optical emission would originate from regions close to the 
pulsar, IC scattering close to might be reasonably approxi- 
mated by the quasi-isotropic expression P c . This would result 
in a possible TeV output Lie ~ riQjMAVP c ~ 10 35 ^ sin a cos a 
erg/s, exceeding the existing upper limits noted above. On the 
other hand, if the NIR-optical emission is produced close to 
7l as proposed in some models (e.g., IPacini & Salvat 1983; 
Crusius-Watz el et al. 20011) . anisotropic inverse Compton scat- 
tering may severely reduce this power output (by up to a factor 
if/ 4 where <A <s 1 is the pitch angle, cf. |MorinTl 981 ), suggest- 
ing a TeV contribution well below the one produced by thermal 
IC. The empirical scaling of the near infrared-optical (and y-ray) 
flux with the magnetic field at the light cylinder seems in fact to 
reinforce a scenario where the emitting region is located close to 
the light surface dShearer & Golden 20011) . 

3.2. Curvature radiation 

Approaching the light surface, field line bending may no longer 
be negligible so that a particle may efficiently lose energy due 
to curvature radiation. In analogy to synchrotron radiation, cur- 
vature radiation can be described as emission from relativis- 
tic charged particles moving around the arc of a circle, chosen 
such that the actual acceleration corresponds to the centripetal 
one (e.g., Och elkov & Usov 19 80). The critical frequency where 
most of the radiation is emitted is given by 



3c 
4jtR c 



r 



(13) 



which for, e.g., R c ~ r L ~ 10 8 cm (cf. Gold [19681 [1969b yields 
v c =! 5 x 10 22 (yl 10 7 ) 3 Hz or a curvature photon energy of about 
0.2 (y/10 7 ) 3 GeV The energy loss rate or total power radiated 
away by a single particle is 



c 3Ry 



(14) 



The characteristic cooling timescale t c = yniQC 1 jP c thus be- 



/, - \mRl[— \\ 

m e I y 3 



(15) 



To find the maximum electron Lorentz factor attainable in the 
presence of curvature radiation, we can again balance f acc (Eq.[3]l 
with t c to obtain 



y max -1.2x 10 9 P 2/5 I — I — 



4/5 



1 



1/5 



• 2/5 

sin ' a , 



(16) 



indicating that for a Crab-type pulsar (P = 0.033 s, R c ~ r\) 
with,. e.g., yo ~ 10 4 curvature radiation constrains achievable 
electron Lorentz factors to y ^ 4.8 x 10 7 sin 2 ^ 5 a. This con- 
firms that for Crab-type parameters, electron Lorentz factors up 



to ~ 10 7 might in fact be obtained. We can again estimate the 
possible curvature output L^ eV at GeV energies, by multiplying 
the particle number with the single particle curvature power P c 
to obtain 



r GeV 



— "GM^f &r*s7 - 



S/s 



(17) 



using AV (with X — 1) as defined above. Note that curva- 
ture radiation in principle results in an additional IC (Klein- 
Nishina) photon target field close to with energy density 
«ph,c ~ L^ eV /(47rr 2 c). For Crab-type values, the resultant en- 
ergy density would be a factor of a few higher than the thermal 
one. Yet, due to both, the further reduced Klein-Nishina cross- 
section and anisotropic scattering conditions, this curvature pho- 
ton field is not expected to significantly modify our IC consider- 
ations above. 



3.3. Pair creation and yy-absorption 

In the magnetosphere of a pulsar, pair creation via magnetic 
photon absorption (y + B — > e + + e~ + B) becomes kine- 
matically possible as the magnetic field can absorb momen- 
tum and thereby ensure momentum conservation. For a pho- 
ton, traveling with pitch angle i/r to the local magnetic field, 
the threshold conditions for this being possible is e p h sin t// > 
2m e c 2 . The absorption coefficient for this process is k =t 1.4 X 
10- 13 (Bsin^/B e ) T{X) [1/cm] where B c = m e c 3 /(eh) =* 4.4 x 
10° G, A = 1.5 (B sin il//B c )(hv/m e c 2 ) and T(A) is the Erber 
function (lErber 1966MTsai & Erber 19741 1 . The function T(A) is 
very sensitive to A. While A might be high near to the stellar 
surface, only moderate values are expected close to the light sur- 
face. For small A <sc 1 one has T(A) 0.46 exp(-4/A) (cf. also 
Daugherty & Hardi ng" 1983[ l, while for large A » 1 one finds 
T(A) ^ 0.9 The location of the absorbing surface may 

thus be approximated by at 1, which gives 



A-' -0.25 ln[6.4 x 10" I4 (B/B C ) r] . 



(18) 



For a Crab-type pulsar (B ~ 10 6 sin 3 o' G close to r\) one 
finds A~ l =i 7.6, so that a photon with energy above the cut-off 
e p h = hv ~ 2/(sin 3 a simp) TeV will undergo magnetic absorp- 
tion. As the pitch angle is usually very small (\J/ ~ 10~ 2 ), we do 
not expect this of significance for Crab-type pulsars at photon 
energies below 50 TeV. 

Apart from one-photon pair production, energetic pho- 
tons may also undergo photon-photon interactions (y + y s — > 
e + + e~) with background soft photons of energy e s (e.g., 
Chiang & Romani 1994). Let us thus consider the following 
cases: 

(1) In the case of TeV photons, the threshold condition requires 
the presence of soft photons with energies e s = 2(m e c 2 ) 2 /e = 
0.5(1 TeV / e) eV or larger. The cross-section for yy-pair produc- 
tion has a sharp maximum of ~ 0.20-7- at 2e s , so that the optical 
depth can be approximated by Tie) ~ Q.\cr T L{2e s ) Ar / (4nr 2 ce s ) 
where L(2e s ) is the corresponding luminosity at which the peak 
occurs and Ar (r - r\) the path length. For a characteristic 
(observed pulsed Crab) photon field of L(O.leV) <K L(leV) ~ 
(1 -3)x 10 33 erg/s (e.g., |Eikenberry et al. 1997llSollerman 20031 
this would result in 



r(e) - 1 



L(l eV[l TeV/e]) 
3 x 10 33 erg/s 



) ( 1 TeV ) 



sino- . 



(19) 
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Fig. 1. Maximum Lorentz factors versus the inclination angle: 
Ymax (fr° m co-rotation) and y^ ax (curvature radiation) are repre- 
sented by the dashed and the solid line respectively. The set of 
parameters emplyed is: M = 5 x 10 3 , P - 0.033s, y = 10 4 . 



noting that (r L /r)(l - r L /r) < 1/4 for r > r L . Equation $19[ thus 
suggests that a substantial fraction of TeV photons (if not all) 
may in fact be able to escape absorption. 

(2) The situation could be somewhat different for the GeV cur- 
vature photons. In this case yy-absorption requires the presence 
of soft photons in the X-ray regime with e s = 0.5 (1 GeV/e) 
keV or larger. In the case of, e.g., the Crab pulsar, the X- 

10 36 erg/s ([Kuiper et al.200T 



ray luminosity is of order L x 

Poss enti et al. 2 002, M assaro et al. 2006t , so that the optica 
depth becomes 



r(e) a 0.34 



L(l keV[l GeV/e]) 
10 36 erg/s 



)(l Gev) 



sin a . 



(20) 



To first order, the soft X-ray flux of the Crab follows 
L(E) a 1Q 36 (10 keV/E)" - 4 erg/s (e.g., IMassaro et al. 20001 
Kuip er et al. 200 1) . Hence, for photons of energy e ~ 20 GeV 
or 50 GeV, the relevant target luminosity would be L(2e s ) - 
1.2 x 10 35 and 8.3 X 10 34 erg/s respectively, implying an opti- 
cal depth r =! 0.8 sin a and r ^ 1 .4 sin a. As curvature radiation 
typically results in a quasi-exponential GeV tail, no significant 
super-exponential suppression is expected to occur at energies 
below several tens of GeV, which seems consistent with recent 
results based on the detection of pulsed emission above 25 GeV 
from the Crab dAliu et al. 20081 1. 



4. Conclusions 

Based on our idealized model analysis above, the efficiency of 
particle acceleration close to the light surface is essentially lim- 
ited by curvature radiation reaction or the Alfven corotation con- 
dition. 

We have studied the implications for conditions applicable to 
Crab-type pulsars, assuming a plasma-rich environment with 
pair density exceeding the primary Goldreich-Julian (close to 
the star) one by M ~ 10 3 . In this case, electron Lorentz factors 
up to y ~ 2 x 10 7 appear possible. Synchro-curvature radiation 
could then lead to a relatively strong (averaged) power output of 
~ (10 34 - 10 ) erg/s at ~ 2 GeV that would be consistent with, 



e.g., EGRET observations of the Crab (Kuiper etal. 2001). The 
emissivity of curvature radiation can be described by 



1 [3 e 2 v f°° e~' 

2 V n R c v c J v/Vc yjx 



(21) 



For a power law-type distribution of particles, for example, 
N(y) = Noy~ s , one can show that for high frequencies (v » v c ) 
the intensity I v behaves as 



fpt. 



y)N(y)dy oc 



(22) 



so that in the case of the Crab the curvature output might be 
expected to decay quasi-exponentially for energies above ~ 10 
GeV, as indeed suggested by recent very high energy gamma- 
ray measurements (Alb ert et al. 20081 lAliu et al. 2008b . Inverse 
Compton (Klein-Nishina) up-scattering of thermal soft photons, 
on the other hand, could result in a (quasi-isotropic) power out- 
put of ~ 10 31 erg/s at TeV energies, consistent with, e.g., ex- 
isting IACT constraints on the observed (non-steady) emission 
from the Crab. 

For older pulsars (e.g., P ~ 1 s) and M » 1, co-rotation usu- 
ally imposes the strongest constraint, so that achievable max- 
imum Lorentz factors are typically limited to ~ (10 4 - 10 5 ). 
Although curvature radiation may then peak in the optical-UV 
(up to ~ 10 eV) and inverse Compton (Thomson) scattering of 
curvature or thermal photons could result in very high energy 
emission up to ~ 50 GeV, their associated power is negligible 
due to the small electron Lorentz factors and the substantially re- 
duced target photon energy density. Hence, within the approach 
considered only young pulsars (P ~ 0.1 s) might be expected to 
produce detectable high energy gamma-ray emission. 

The proposed scenario could in principle work for a variety 
of angles, so that the resultant emission needs not necessarily to 
be strongly pulsed. Note that for most circumstances, the major 
condition limiting the Lorentz factors of electrons results from 
co-rotation. Yet, for large inclination angles, curvature radiation 
reaction can become dominant over co-rotation, cf. Eqs. (0]i and 
([Tol l and see Fig.Q]for illustration. 

The analysis presented is based on a number of idealizations, 
which we plan to remedy in future studies. This particularly in- 
volves the assumptions of, e.g., quasi-straight field lines and a 
single particle approach in which plasma effects are neglected. 
On the other hand, one of the strengths of the present concept 
is its ability to explicitly take inertial effects into account and so 
to allow to estimate the size and extent of the VHE regions in 
young pulsars. 
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